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ABSTRACT 

We study the (i-dimensional random Ising model using a suitable type of Bethe- 
Peierls approximation in the framework of the replica method. We take into account 
the correct interaction only inside replicated clusters of spins. Our ansatz is that the 
interaction of the borders of the clusters with the external world can be described 
via an effective interaction among replicas. The Bethe-Peierls model can be mapped 
into a single Ising model with a random gaussian field, whose strength (related to the 
effective coupling between two replicas) is determined via a self-consistency equation. 
This allows us to obtain analytic estimates of the internal energy and of the critical 
temperature in d dimensions. 

PACS NUMBERS: 05.50.+q, 02.50.+S 
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Introduction 

The mean field solution and its improvements, such as the Bethe-Peierls 
approximation [1,2], give good approximations of the critical temperature and of 
the internal energy in many statistical models. We show that the same methods 
can be applied to spin glasses by considering the overlap among replicas instead of 
the magnetization. Indeed, the appropriate ansatz for spin glasses is assuming that 
the effect of the thermal bath on a replicated cluster of neighbors spins produces an 
effective coupling among replicas. We shall give an a posteriori justification of such 
a hypothesis by proving that the Sherrington-Kirkpatrick (SK) model [3] is recovered 
in the limit of infinite dimension. 

This paper considers the Ising model with independent random nearest neighbor 
coupling Jij in absence of external magnetic field. Our main result is that, in the 
Bethe-Peierls approximation, this model is equivalent to a single Ising model with 
a random gaussian field whose strength is related to the effective coupling between 
two replicas. We thus obtain an estimate of the internal energy and of the critical 
temperature in any dimension. 

In section 1, we introduce the Bethe-Peierls ansatz for spin glasses in the framework 
of the replica method. 

In section 2, we prove that this ansatz leads to the SK model when d — > oo 

In section 3, we show that, under the hypothesis of no replica symmetry breaking, 
the d-dimensional model can be mapped into a single Ising model with random 
gaussian field. We also explicitly compute the replica symmetry solution for the 
internal energy in d dimensions for the spin glass with dichotomic random coupling 
J = ±1. 

In section 4, we show that our method allows to compute in a simple way the 
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critical temperature Tc{d). The result is very accurate at high dimension. 

In section 5, we discuss the possibility to use our ideas to implement a clever 
numerical scheme for determining internal energy and critical temperature of d- 
dimensional spin glasses. 

1. Bethe-Peierls ansatz for spin glasses 

The partition function of the Ising models on a lattice of N sites with nearest 
neighbor couplings which are independent identically distributed random variables, 
in absence of external magnetic field, is 

Zn{P, Uij}) = ^ n eMPJij<^i<^j) (1-1) 

{s} (iJ) 

where the sum runs over the 2^ spin configurations {s}, and the product over the dN 
nearest neighbor sites 

In the thermodynamic limit almost all disorder realizations has the same free 
energy, i.e the quenched free energy 

1 



where A indicates the average of an observable A over the distribution of the random 
coupling P{Ji,j) . In the following we assume that the P{Jij) is such that Jij = and 

'Jij - -L- 

On the other hand, it is trivial to compute the so-called annealed free energy 

^" = -JL- i^'"^. (1-3) 

corresponding to the free energy of a system where the random coupling are not 
quenched but can thermalize with a relaxation time comparable to that one of the 



spin variables. For instance, in the case of dichotomic random couphng = ±1 with 
equal probability, one has 

/a = (ln2 + dlncosh/3) (1.4a) 

while for gaussian coupling, i.e. P{Jij) = C^*^/^ / , one has 

/« = -/?-^(ln2+ ^) (1.46) 

However, is in general very different from the quenched free energy. In order 
to compute (1.1), it is convenient to use the replica trick [4]. Let us thus consider 
n non- interacting replicas of the disordered system labelled by a = l,---,n. The 
corresponding partition function is 

Z- = J2e.p\pj:J2JiM-^a^A (1.5) 

{s} \ a=l(i,j) J 

where the sum runs over the 2^^ spin configurations {s} of the replicas, 

{s} = {.W},..., {.(-)} 

with 

After having performed the average and found an analytic continuation at real 
n-values, the quenched free energy is given by 

In Z = lim i In (1.6) 

n^O n 

Even in two dimensions, there is no exact solution for this problem. The first 
non-trivial approximations of the quenched free energy can be obtained either by 
constrained annealed average [5] or by improved mean field approximations of the 
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Bethe-Peierls type. Recently we have introduced such an approximation in the dual 
lattice made of square plaquettes in two dimensions [6]. However, there is no solution 
of the self-consistency equation at low temperature and it is not trivial to generalize 
the approach at higher dimensions. For systems with diluted quenched disorder, a 
different type of improved Bethe-Peierls approximation (the Cluster variation method) 
has been studied in [7] without using the replica approach. 

In this paper, we want to work directly on the real lattice, by taking into account 
the correct interactions inside a pile of replicated cluster made of a central spin ctq 
and of its 2d nearest neighbor {ufc}, and by considering only an effective interaction 
with the external world. Note that in the 2d case, the clusters are crosses made of 5 
spins. 

Separating the two contribution (crosses plus external world) in the partition 
function, we get 

(n 2d n 
a=lk=l {sext} a=l (ii)#(Ofc) 

(1.7) 

where 

Jk = Jok 

are the coupling between the central spin of the cross and its neighbors on the border. 
The first sum in (1.7) runs over the 2^'^'^'^^^'^ spin configurations {scr} of the replicated 
crosses labelled by (uq , erf , • • • , 0-2^ ) with a= 1, . . . , n while the second sum over all 
the other spins. The expression obtained by computing the second sum depends 
only on the 2dn lateral spins (j^. The correct Bethe-Peierls ansatz is given by the 
assumption that the interaction among the lateral spins of the replicated crosses 
and the external word forces an effective interaction among different replicas with 
a constant that should be determined via a self-consistency equation. In other 
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terms, our Bethe-Peierls ansatz is 



exp [/3 X: E J^o^^ I = ^(/^) f E ^"/^ E I (1-8) 

{Sext} V a=l(ij^Ofc) / \a>/3 i=fc / 

where K{j3) is a multiplicative constant which depends on the temperature but not 
on the lateral spins. One expects that //q:/? = in the high temperature phase, while 
it must have a non-zero value in the glassy phase. 

Therefore, instead of the (1.7), we have to compute an effective partition function 



Z. 



n 

2d 



= E n ^-p E J'^^o''^ h-p E /^^/^^"^ (1-9) 

{sc,.}fc=l \ a=l / \a>/3 J 

A further simplification can be reached for dichotomic coupling = ±1 where 
one can perform the gauge transformation — > Jk cr^^ on the lateral spins, leaving 
the free energy unchanged. In this case the averaged partition function (1.9) becomes 

= E n ^-P (/^E^J"^^"^) exp (j^f^aA^'^aiA (1.10) 

This relation implies the rather surprising result that a non-disordered Ising system 
exhibits the same behavior of a spin glass if one imposes the appropriate interaction 
among different replicas. 

At this point, the effective coupling fJ''^i3{P) is given by the self-consistency equation 

<^P)n = lim (1.11) 

where represents the thermal average over the replicated system. Then, the 
Bethe-Peierls estimate of the internal energy is 



UBpiP)= lim-;l 



d 



(1.12) 



Let us anticipate that the Bethe-Peierls approximation predicts a phase transition at 
a critical temperature Tc{d) above which //*^ = 0. As a consequence the Bethe-Peierls 
solution coincides with the annealed one in the high temperature phase, i.e. 

Ubp = ^ [PfaiP)] for p<Pc (1.13) 
2. The infinite dimensional limit 

The model defined by (1.9) becomes the infinite range SK model [3] in the limit d — > oo. 
This result has a great importance since provides a good evidence that we have chosen 
the correct Bethe-Peierls ansatz for spin glasses. In this section, we prove that the 
self-consistency equation (1.11) in the limit (i — > oo gives the equation for the overlap 
of the SK model. 

Let us recall that the averaged partition function of the infinite range SK model 
after some simple algebraic manipulation becomes 

N 



{zy 



{a} \oi>l3 



(2.1) 



where the sum is on the 2"^ realizations {a} of the n spins a^^\ . . . ,a^'^\ In the high 
temperature phase T > Tc, one has = so that Z"' = (Z)"^, while in the glassy 
phase one has a non trivial overlap qctp = l^piT) which maximizes {Z)'^^ that is 

QaP = (aV^) = y ^ r ^ (2.2) 

In order to get the correct d — > oo limit of the self-consistency equations (1.11), 
we should use the rescaling 

P -J= IJ^aP P^l^ap (2.3) 

y 2a 



Now, the disorder average in Z„ is easily performed since at large d the first exponential 
in (1.9) can be expanded in Taylor series up to the second order, so that 



exp(/3 {2d)-y^JkSk = 1 + /3(2d)-V2 j^^fe + ^ Si + 0{d-^/^) = 



= 1 + ^Sl + 0{d-'/') = eM^/l) + 0{d-''') 

where = X^a '^o'^^ ^^if^ • The distribution of the coupling is irrelevant provided that 
J = and = 1. Therefore, a part small correction 0((i~^/^), the partition function 
(1.9) becomes 

2d 2d \ 

E -o"' <^i« Yd E •'i"' "i" + E ' -i"' (2-4) 

a>^ fe=i y 

implying that 

V ^(") ^(^) f^^T. T/^s „ ^(a) 

< 4"' >= ^ ^ ;>-p('^£»>/>^^ '-^ j ^ ^) 

E{afc} exp Ea>/3 /^a/3(^r 4 j 

where /c is one of the lateral sites of the d-dimensional crosses and the sum is on the 
2" realizations crfe of the n spins cr^^\ ■ ■ ■ , c^"^. 

On the other hand, in the limit (i — > oo, the corresponding relation for the central 
spins can be written as 

since one has 

^ 2d 

/ (a) (13) ^ r (a) (/9) /o 

k=l 

A direct comparison of (2.6) and (2.5), shows that the self-consistency equation (1.11) 
is satisfied only if 



/ (a) (13) ^ 



that is the equation for the overlap of the SK model. We can thus identify the coupling 
Hais with the overlap qa/3 for (i — > oo. 

3. Replica symmetry solution in the Bethe-Peierls approximation 

It is possible to obtain the replica symmetry solution of a d-dimensional spin glass 
in the Bethe-Peierls approximation. We must note that in the case fia/3 = A*) the 
averaged partition function (1.9) of n replicated crosses is 



2d / n 

- E n -P E -i"^ ) exp I ^ E 1 (3-1) 

{sc,.}A;=l \ a=l / y a>/3 

A part constant multiplicative factors, it can also be written as 

2d 

■exn(/9.7.V. exn^ 



= E n ^-p E -o'^^ -t' I ^-p f I E -t' ) (3.2) 



that is bilinear in (T]~ . In order to linearize (3.2), we should use the standard gaussian 
identity 

exp(a;^/2) = / dw exp(— a;^/2) exp(a;a;) 

v27r Joo 



so that one has 

2d 



= E n ^^p E( ^0 + vi^^k 4"^ ) (3-3) 

{Scr} k=l \q:=1 / 

where 

-, „ „ 2d 2d 

^= (2^ J '" J Udi^ke-'^^'^'IlPWdJki^ (3.4) 

indicates now the average over the standard gaussian variables and over the 
coupling Jfc between central and lateral spins. 

This transformation has the advantage to allow for a factorization of the product 
over the replicas in (3.3), implying that 

(3.5) 



with 

2d 

$ = JJ exp (/3 Jfc (Jo cTfc + ^WfeCTfe) (3.6) 

This is the main result of the section. It estabhshes that, in the Bethe-Peierls 
approximation, the rephca symmetry solution is equivalent to that one of a single 
Ising model with a random gaussian field applied to the boundaries of the cross. This 
field has a strength related to the coupling among replicas and describes the interaction 
of the cluster of 2(i + 1 spins with the external world. 
The explicit sum over the lateral spins (jfc's gives 



Zn = \Y. ^m(^o) (3.7) 

\{Scr} / 

where is the non-normalized weight of the central spin 

2d 

Wt,{ao) = 2cosh(/?Jfc(7o + ^oJk) (3.8) 
fc=i 

obtained after summing over the configurations of the 2d lateral spins Ufe. The 
probability of the central spin thus is 

""'^^'^ ^ W.iao = ^fwliao = -1) ^'"'^ 
and is itself a random quantity depending on the 2d random gaussian fields and to 
the 2d random coupling J^. 

Because of the replica symmetry, the self-consistency equation (1.11) for 
determining //*, and so the needed strength of the random field, assumes the simpler 
form 



< (To >2 = < (Ji >2 (3.9) 

where the thermal average of the central spin is 

<ao>= Yl ^oJ"/.(cTo) (3.10) 

o-o=±l 
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and the thermal average of one of the lateral spins is 

<(Ti>= ^ tanh(/?JfcC7o + ^/^luJ^) P^,{(Tq) (3.11) 

ao=±l 

In order to find the internal energy we have to compute 

^ 2d 

lim — InZn = In IT 2 cosh(/3 uo + J]lijOk) (3-12) 

(T0=±1 fc = l 

and then, following (1.12), the internal energy is given by a derivative at = //*, 
solution of (3.9) 

UspiP) = -d JkCTQ tanh(/3Jfecro + ^/JIuJi) -P/i(<7o) (3.13) 

It is worth stressing that in the limit d ^ oo, the self-consistency equation (3.9) 
becomes, by virtue of the results of section 2, 

11= lim < (Ti >2 (3.14) 

d— >CXD 

The above expression, after performing the rescaling (2.3) gives the replica symmetry 
solution for the overlap of the SK model in the glassy phase, 

/x = tanh^(/3a;V^) (3.15) 

where uj is again a standard gaussian. 

For the ±J model, after the gauge transformation Jkcrk crfc, the probability 
-F/i(co) depends only on the gaussian fields and is independent of the coupling J^. We 
can thus put Jfe = 1 in the formulas from (3.3) to (3.13) and the average (3.3) should 
be taken only over the 2d gaussian variables a^^'s. 

Fig 1 shows the replica symmetry solution T'^/j,* (T) /2d as function of the rescaled 
temperature T/ \/2d at d = 2 , 3 , 4, 6 for the ± J model. The effective replica coupling 
11* vanishes above the critical temperature Tc{d) as we shall discuss in the next section. 
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As a consequence the internal energy Ubp is equal to the annealed internal energy at 
T>Tc. 

Another peculiar effect arises in the low temperature regime where //* (T) 
decreases below Tq pa 0.5. This vanishing of the coupling IJ'*{T) can be understood 
by the following qualitative arguments. Consider the model defined by (3.6) with Jj. = 
1 which is originated by the ±J model. On the lateral spins, there is a competition 
between the coupling with the random field ^/JiukO'k and the ferromagnetic interaction 
/3(Tocrfe. It gives origin to a frustration of the system below Tg. However, if the 
temperature is very low, the ferromagnetic interaction dominates since the work 
/ (3 necessary to win the tendency of the spin ak to align with the field vanishes. In 
correspondence the system would become ferromagnetic with a ground state Uq = —d 
as it happens for the annealed model. Such a regime is clearly unphysical, and one can 
trust in our results only when the work made to destroy the long range order 

in the glassy phase is a non-increasing function of the temperature, i.e. for T >Tq. 

Notice that the mean field coupling is obtained at increasing d, by the 'a posteriori ' 
rescaling fx* ^ fx* / (2d/?^) for the coupling and T ^ T' = T / \pld for the temperature 
in the sequence of lines for the d-dimensional functions. It is therefore natural our 
choice for the coordinates of Fig 1. Let us also stress that the rescaled temperature 
below which T^//* decreases is Tg = 0.5 / ^J2d so that the unphysical ferromagnetic 
regime disappears when cZ — > oo. 

The internal energy \J b p (T) /dis shown in Fig 2 for the ± J model at d = 2 , 3, 4 , 6. 
An estimate of the ground state energy Uq can be obtained by Ubp{Tq) as previously 
argued. Using this hypothesis we get 

Uo = -1.51 at d = 2, 
Uq = -1.88 at d = 3, 
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Uq = -2.204 at d = 4, 
Uo = -2.718 at d = 6 

At d = 2 we can compare our analytic estimate with the numerical result [8], 
Uq = -1.404. 

It is an open issue to understand whether better estimates can be obtained via 
(1.12) with a replica symmetry breaking solution A*a/3- 

4. Phase transition and critical temperature in finite dimension 

The Bethe-Peierls method and its improvements are able to give accurate estimates of 
the critical temperature in a disordered system. Using a replica symmetry approach, 
good analytic results have been obtained for diluted spin glasses [9] and other randomly 
frustrated systems with finite connectivity [10]. 

In the framework of the results of the previous section, we should note that at the 
transition point, the order parameter ji* vanishes. Therefore, the critical temperature 
can be computed from (3.3) considering only the first order of its expansion in jj.* . 

Let us first compute the thermal average of the central spin 



Since this expression appears in (3.3) only in a squared form it is not necessary to 



2d 

< (70 > = -\/a* tanh(/3 J/.) ^a;/- + 0{ii) 

k=l 



(4.1) 



compute higher orders than y^/I. Analogously, the thermal average of one of the lateral 



spins IS 



2d 

< (Ji > = tanh(/9 Ji) ^ tanh^(/5 Jj.) cuk + \^ + 0(//) 

fc=2 



(4.2) 



Inserting this expression in the consistency equation (3.9) one obtains: 




(4.3) 
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where J is one of the couphngs. Eq. (4.3) gives the critical temperature T(. = (3^^ as 
a function of the dimension 

tanh^(/3eJfe) = (4.4) 

In the case of the ± J model, this equation becomes 

tanh2(/3e) = (4.5) 

In Fig 3, we compare the Bethe-Peierls critical temperature (4.5) with the numerical 
result obtained in the literature for the ±J model [11]. 

We have also computed the critical temperature of the gaussian model via a 
numerical solution of (4.4). In this case: 

Tc = 1.19 at d= 2, 

Tc = 1.81 at d = 3 (numerical result = 1.0), 
Tc = 2.28 at d = 4 (numerical result Tc = 1.8), 
Tc = 2.67 at d = 5, 
Tc = 3.06 at d = 6. 

Let us remark that Tc is finite in d = 2. This spurious transition is a typical 
and well-known effect of mean field approximations. In fact. The Bethe-Pcicrls 
approximation, gives a lower critical dimensionality = 1, where Tg = 0, while 
there is a good numerical evidence that dc = 2. On the other hand, the higher the 
dimensionality, the better our estimates. In the limit of infinite dimension, after the 
usual rescaling (2.3) of the temperature, from (4.5) one obtains /?c = 1 which is the 
critical temperature of the SK model. 

It is possible to improve the estimate of the critical temperature in a systematic 
way by considering larger cluster instead of a cross made of a single central spin ctq 
and of its 2d neighbors cr/c, as we shell discuss in the conclusions. 
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For instance, we have considered a plaquette of four spins plus the 8{d—l) spins 
that are their nearest neighbors in the ± J model. Applying our Bethe-Peierls ansatz 
to the replicated plaquettes, the equation for the critical temperature Tc{d) is again 
given by the solution of a rational function oft = tanh(/3). After a lengthy but trivial 
calculation one has 



l + {2d-3)r -2{d-l)r + 2{d-l){r -r) t 



(1+^2)2(1+^4)2 (1+^4)2 



= 
(4.6) 

In this case, the lower critical dimension is dc = 15/11 instead of cZc = 1 found for the 
cross. The critical temperature Tc{d) obtained by (4.6) is shown in fig 3, too. 

In our opinion, looking at increasingly larger clusters it is possible to determine the 
critical temperature of a d-dimensional spin glass as the zeros of rational functions of 
t = tanh(/?) reaching an accuracy much larger than that one given by direct numerical 
methods. Moreover, one can also hope to find a converging sequence of lower critical 
dimensions, simply considering the zeros of the rational functions with t = 1 (i.e. 
Tc = 0). 

5. Conclusions and perspectives 

The properties of finite dimensional Ising spin glasses are largely unknown. The lower 
critical dimension itself is not known although most of numerical simulations indicate 
d — 3 as the lowest dimension which exhibits a glassy phase at finite temperature. 
Furthermore, even if the glassy phase is present, the existence of replica symmetry 
breaking at low dimensionality is still controversial. All that is a clear indication of 
the difficulties encountered when one tries to extract information directly from the 
model. 

In our approach we simplify the task. Indeed, when we assume replica symmetry, 
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our approximation reduces itself to the study of the model 

2d 

^= n exp (/3 Ja; cJo (Tfc + ^Jjiukcrk) (5.1) 

where both the ujk and the are quenched random variables. It should be noticed 
that in this model, one only deals with 2d+ 1 spins and Ad quenched variables at most. 
The model is completed by the self-consistency equation (3.9) that we rewrite here as 



<cTo>'= ; (5.2) 

2^k 

where the second term of (5.2) is the mean of the overlap on the lateral spin. 

The validity of our approach stems from the possibility of a systematical 
improvement. Following a standard technique we can replace (5.1) by 

Z = ^Ylexp{/3 Ju> ai ) exp {/3 Jik o-^ cTfc) exp {^ujkcrk) (5.3) 

{s} 1,1' l,k k 

where the first product is on all the first neighbor spins of a hypercube, the second 
product is on the couples of spins formed by lateral spins labelled by k on the faces 
of the hypercube and their first neighbors, the third product is simply on the lateral 
spins. This model is completed by the self-consistency equation 



where the first sum runs on all the spin of the hypercube and the second one on the 
lateral spins. The linear dimension of the hypercube can be progressively increased, 
and one expects to converge to the right result in the limit of large hypercubes. In our 
opinion, this might be a powerful numerical tool to determine the internal energy and 
the critical temperature of a spin glass, superior to a direct approach by Montecarlo 
simulations. 

16 



Let us also mention the major open problem from a theoretical point of view. It 
is the search of the solution of the Bethe-Peierls equations with replica symmetry 
breaking, to see whether the unphysical behavior of the internal energy at low 
temperature disappears as it happens in the Parisi solution [12] of the SK model. 
A first step can be reached by looking for a solution with only one breaking. This can 
pave the way to the comprehension of the glassy transition at finite dimension. 
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Figure Captions 



Fig. 1 Replica symmetry solution of the self-consistency equation fi*/{2d/3'^) as function 
of the rescaled temperature T/ \f2d for the ± J model at (i = 2, 3, 4, 6. The larger 
the dimension, the higher the corresponding line. The dashed line indicates the 
infinite dimensional limit (overlap of the SK model). 

Fig. 2 Annealed internal energy Uajd = tanh(/3) (dashed line) and the Bethe-Peierls 
solutions Usp/d (full lines) versus temperature T = (3~^ for the ±J model at 
d = 2, 3, 4, 6. The larger the dimension, the higher the corresponding line. The 
dotted lines are the estimates of the ground state energy obtained by imposing 
that jjL* is a non-decreasing function of the temperature. 

Fig. 3 Rescaled critical temperature Tc/\/2d versus the dimension d. The Bethe-Peierls 
solution for the ± J model given by (4.5) is indicated by a full line. The improved 
estimate obtained by (4.6) where the cluster is a plaquette of four spins instead of 
a central spin is indicated by a dashed line. The squares are the numerical values 
of Tc/ \/2d for d = 2, 3, 4, 6 joined by a dotted line. 
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